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Abstract 

We use the 1/N expansion formalism in a systematic study of high-spin 
states in the sd and sdg boson models with emphasis on spin dependence of 
moment of inertia and E2 transitions. The results are applied to the high-spin 
states in the actinide nuclei 232 Th, 234 - 238 TJ ) w h ere the need for g bosons is 
especially acute but until now, no realistic calculation existed. We find that 
the (i-boson energy plays a crucial role in description of the high-spin data. 
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One of the early predictions of the the interacting boson model JT|] (IBM) is the 
so called boson cutoff effect. Due to finite number of bosons N in the system, the 
ground band terminates at spin L = l ma xN, where l max denotes the maximum boson 
spin. As a result, yrast B(E2) values are predicted to fall off relative to the rotor 
values, vanishing at the maximum spin. In the sci-IBM, the cutoff occurs at L = 2N 
which is low enough to be experimentally accessible and, therefore, has prompted 
many investigations (see, for example, [[| [3J [|, |^, 0] ) . When the predicted falloffs did 
not materialize in any of the experiments this was interpreted as due to insufficient 
collectivity of the sd-boson system at high-spins which could be ameliorated by 
including g bosons 0, ||, @. Since then, evidence for g bosons from low-spin spectra 
have continuosly grown with concurrent developments in the sdg-IBM calculations 
(see [|K], [Tl| for recent reviews). Ironically, the high-spin data, which have been the 
harbinger of g bosons, could not be treated in realistic sdg-IBM calculations. This 
happens because the basis space for deformed nuclei (N = 12 — 16) is too large and 
has to be severely truncated for diagonalization, limiting their validity to low-lying 
states. Although the SU(3) limit has been used in discussing some features of high- 
spin states in the sdg-IBM ^ [|, realistic calculations indicate that this limit 
is rather strongly broken Hence it can only provide a qualitative picture but 
otherwise is not very useful in applications to spectroscopy, or in addressing some 
basic questions on the shortcomings of the s<i-IBM raised by Bohr and Mottelson 
13 ], e.g. the spin dependence of moment of inertia. An unfortunate consequence of 



the lack of sdg-IBM calculations is that the sd-IBM results are still being used in 
comparisons with high-spin data with obvious negative connotations (see Refs. || 



and |14]] for recent examples in experimental and theoretical studies). 

The angular momentum projected mean field theory, which leads to a 1/N ex- 



pansion for physical quantities |15] , presumably offers the only viable alternative for 
a realistic description of high-spin states in the sdg-IBM. An obstacle in realizing 
this goal, namely evaluation of matrix elements up to order (L/N) e which is neces- 
sary for an accurate representation of high-spin states, has been recently overcome 
through the use of computer algebra |LJ|. The purpose of this letter is to use the 
extended 1/N expansion formulas in a systematic study of high-spin states (before 
the backbend) in boson models. We address, in particular, the moment of inertia 
and boson cutoff problems in the sd-IBM, and discuss how far they are resolved 
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with the inclusion of g bosons. Applications are made to the actinide nuclei, 232 Th, 
234_238 U, for which the ground band has been followed up to spins L > 2N, and 
hence the need for g bosons is most acute. 



The 1/N expansion method has previously been discussed in detail [15| and the 
recent extensions to higher orders are given in Ref. [HJ. Therefore, we give only a 
short account of the formalism here. The starting point of the 1/N calculations is 
the boson condensate 



|iV,x) = (iV!)- 1 / 2 (6t)^|0), &t = £ X/6 



t 



where b l0 denote the boson creation operators and X\ the associated mean fields. 
The summation index runs over I = 0, 2 in the sd model and I = 0, 2, 4 in the sdg 
model. The axial symmetry assumption in Eq. (1) has been justified in Ref. [16[ 



by comparing the 1 /N expansion results with the exact diagonalization ones. For a 
given Hamiltonian H, one evaluates the projected energies El = (N,x\HPq \N,x) 
and determines x\ by variation after projection (VAP) |T5). The resulting energy 
expression is a double expansion in 1/N and L = L(L + 1), and has the generic form 



^=^ 2 E7^(4^r> (2) 

w {aN) m ^a 2 N 2 ^ 



where a = J2i an d the expansion coefficients e nm involve various quadratic forms 
of the mean fields X\. The coefficients e nm have recently been derived up to the 
order L 3 /N 6 Hl6| . Since the N dependence is not relevant to the discussion, we will 
suppress it for simplicity and rewrite the energy formula @ as 

E L = e x L + e 2 Z 2 + e 3 L 3 , (3) 



where the coefficients e n can be read off from the expressions given in Ref. JTB| . 
We note that similar expressions are used in the geometrical model analysis of 



deformed nuclei [17]. The difference between the two models is that in the IBM the 
coefficients e n follow from an underlying Hamiltonian (which is used in describing 
other properties) whereas in the geometrical model they are directly extracted from 
the data. The moment of inertia problem raised in Ref. |K| refers to the fact that 
in the sd-lBM, i) the e\ coefficient gets a substantial contribution from the dipole 
interaction, L ■ L, which has no dynamical content, and ii) the e 2 coefficient is much 
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smaller than the experimental values. The two problems are in fact interrelated. 
Although the latter can be resolved by renormalizing the moment of inertia at high- 
spins (e.g. by modifying L • L — > L ■ Lf (1 + fL ■ L) flBfl), such modifications are 
purely kinematical in origin and do not address the dynamical problem. 

Another observable of interest in the study of high-spin states is the yrast E2 
transitions which have the generic form 

(L + 2 || T(E2) || L) = aNL(L020\L + 2 0)[mi + m 2 L(L + 3)] (4) 

where a is an effective boson charge, L = [2L + l] 1 / 2 and the coefficients m n are 
given in Ref. [16|. The first term in (||) gives the familiar rigid-rotor result. The 



second term is negative and is responsible for the falloffs predicted in E2 transitions. 

In order to compare the sd and sdg model predictions, we first present a brief 
study of the sd-lBM results. We use the standard Hamiltonian 

H = -kQ ■ Q + k'L ■ L + e d n d , (5) 

where L and n d are the angular momentum and <i-boson number operators, and the 
quadrupole operator is given by 

Q = [s^d + Ss]^ + x[d ] d\ {2) . (6) 

Here brackets denote tensor coupling of the boson operators and bi m = (— l) m 6;_ m . 
For consistency, the same quadrupole operator is used in the E2 transition operator, 
T(E2) = aQ (|4]) as in the Hamiltonian ([5]). Since L ■ L does not play any role in 
the dynamics of the system, we will not discuss it further (it can be easily restored 
by changing e± — > e\ + k'). In presenting systematics, we find it convenient to use 
ratios which eliminate the undesired effects of the scale parameters k and N. The 
energy scale can be fixed, for example, by fitting k to the excitation energy of the 7 
band, E 7 . In Fig. 1, we show four such quantities as a function of q = x/Xsu3 for 
various values of rjd = Ed/Nn. We comment on their behaviour and contrast them 
with the experimental data below. 

a) E-y/Nei: This ratio relates the energy scales of the 7 and ground bands, and 
its mismatch with experiment has been a source of criticism [Q. It is around 4-5 
in the rare-earth region and increases to 8-10 in the actinides. The SU(3) limit 
(q — 1, T)d — 0) is seen to give the maximum value which overestimates it by a factor 



of 2-4. It decreases rapidly with e d and q though, and through a judicious use of 
these parameters, it should be possible to describe this ratio (and hence the the 
moment of inertia) without using the L ■ L term. 

b) N 2 e2/e\\ This ratio measures the deviation from the rigid rotor behaviour (SU(3) 
limit) due to loss of pairing. It ranges from about -0.2 in the rare-earth region to -0.1 
in the actinides. Clearly, it can not be explained by the standard s<i-IBM Hamil- 
tonians currently in use for deformed nuclei which assume r] d = 0, q ~ 0.4 — 0.5. 
However, it is quite sensitive to r\ d values and the experimental range can be easily 
attained by including the <i-boson energy in the Hamiltonian. 

c) N^e^/ei: There is some uncertainty in the extraction of this ratio from data, 
especially in the rare-earth region. In the actinides, it is about 0.01. It depicts even 
more sensitivity to both q and r] d , and therefore its description should not pose any 
problems. 

d) N 2 mijm\. As there is no boson cutoff effect, experimentally this ratio is consis- 
tent with zero. For r] d = 0, it remains rather flat at the SU(3) value which gives 
the maximum possible effect. Introduction of one-body energy, however, reduces it 
substantially, becoming more in line with experiments. 

The (i-boson energy has been mostly neglected in studies of deformed nuclei, 
presumably due to the success of the consistent-Q formalism (variable x with e d = 0) 
in explaining the energy and E2 transition systematics of low-lying states ||12|| . In 
fact, for small values (rjd ~ 1), its effect on low-lying states is negligible and it is not 



really needed in their description |L9j . The above analysis indicates that breaking of 



the SU(3) limit by either the pairing interaction [13 1 or by varying the x parameter 
|T2|j does not lead to a soft enough energy surface which is the main reason for the 
perceived moment of inertia problem in the srf-IBM. The obvious way towards a 
softer energy surface is to include the d-boson energy in the Hamiltonian which is 
seen to vastly improve the description of the spin-dependent terms in the ground 
energies and E2 transitions. 

We next present a similar study in the sdg-lBM. A minimal extension of the 
sd-lBM to the sdg model can be achieved by including the g-boson energy term, 
£ g n g in the Hamiltonian (||), and modifying the quadrupole operator (|6|) to 

Q = [ a tJ + d t § f) + q22 [Sdf) + q2 ,[S~g + g^f) + q4A tfgp). (7) 
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The quadrupole parameters g 2 4 and q^ play an important role in the description 
of the hexadecapole bands but otherwise, the ground band properties are not very 
sensitive to their variations. In order to limit the number of parameters, we scale 
the three quadrupole parameters #22, ?24, 944 from their SU(3) values with a single 
factor q as suggested by microscopies [p0[| . 111 Fig. 2, we show the ratios in Fig. 1 
as a function of q for various values of r] g = e g /Nn with rjd = 0. Before commenting 
on specific ratios, we point out some general features. For large r] g , the g bosons 
decouple and the results converge to those of the sd model shown in Fig. 1. This 
convergence is apparent from the overlap of lines in Figs. 2-a and b but requires 
even larger values of r] g in the case of 2-c and d. It is harder to pin down realistic 
values for r\ g due to lack of data, nevertheless, we quote the literature values for 
comparison which range from 3-6. 

a) E 1 /Ne\: Inclusion of g bosons increases this ratio which is contrary to the ex- 
perimental trend. However, for realistic rj g values, this adverse change is too small 
to worry about. 

b) N 2 e2/ei. This ratio also increases (in absolute value) which is good but again 
too small for realistic r] g values to have any impact. 

c) N A e^/ei: This ratio shows some sensitivity to g bosons, however, it is nowhere 
near the effect of % in Fig. 1-c, and therefore not likely to have much relevance. 

d) A^ 2 m 2 /m 1 : The boson cutoff was the original reason for the introduction of g 
bosons and it is clear from this figure why. In the SU(3) limit, this ratio is reduced 
by a factor of 4 compared to the s<i-IBM. Its q and r] g dependence, however, is 
opposite to that of Fig. 1-d, and things get worse away from the SU(3) limit. For 
realistic parameters, the reduction from the s<i-IBM result (with rja = 0) is less than 
40% which is certainly not enough, and one needs the <i-boson energy to reduce it 
further. 

The somewhat surprising conclusion of the above systematic study is that intro- 
duction of the g bosons, though necessary to describe states with L > 2N, hardly 
improves the dynamics of the boson system. The problems attributed to the sci-IBM 
are, in fact, due to not having a soft enough energy surface and can only be resolved 
by including the ci-boson energy in the Hamiltonian (and not by introduction of g 
bosons alone). 

In the light of the systematic trends discussed above, we carry out fits to the 
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actinide nuclei, 232 Th, 234 - 238 U ; for which the ground band has been followed up 
to spins L > 2N, and hence the need for g bosons is most acute but no realistic 
IBM calculation yet exists. The parameters used in the fits are shown in Table 1. 
The ground band energies (normalized with L(L + 1) for a fairer representation of 
all spins) are shown in Fig. 3. A good description of the ground band energies is 
obtained for all four nuclei with relative errors of about 1%. It is worth emphasizing 
that the moment of inertia and its spin dependence are explained without appealing 
to the L ■ L term or its phenomenological variations. To appreciate the importance 
of the ci-boson energy in the fits, we note from Fig. 2b that with rjd = 0, the slopes 
of the curves in Fig. 3 would be an order of magnitude smaller, leading to almost 
flat lines like in a rigid rotor. In Fig. 4, the calculated E2 transitions along the 
ground band are compared to experiment. Again a reasonably good description of 
the data is obtained. Finally, to show that the ground band results are not achieved 
at the expense of other bands, we compare in Table 2 the predictions for the /3, 7 
band head energies and E2(2p n — > g ) matrix elements with experiment. Note that 
these calculations are the leading order 1/N expansion results and hence are correct 
to order 1/N. The good agreement obtained in Table 2 verifies that the inclusion 
of the (i-boson energy does not detract from the usual quality of the (3 and 7 band 
systematics achieved in the IBM calculations. 

In conclusion, we reiterate that the perceived problems with the sci-IBM in its 
description of spin dependent quantities is not due to lack of higher spin bosons 
but rather due to the energy surface not being soft enough. Inclusion of the d- 
boson energy, together with the extension to sdg space, can successfully resolve 
these problems as demonstrated in the fits to the actinide nuclei. The g bosons are 
necessary for extending the model space but otherwise they play a marginal role in 
the dynamics of the ground band and can not resolve alone the problems mentioned 
above. Details of the present work, including applications to the rare-earth nuclei, 
will be presented in a longer paper |22[ . 
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Table 1: Parameters used in the sdg-IBM calculations, k is in keV and a in eb. 
Nucleus N k q Vd Vg a 

232 Th 12 14.3 0.7 1.85 4.2 0.210 

234 U 13 15.6 0.7 1.62 3.1 0.195 

236 U 14 16.2 0.7 1.62 3.1 0.195 

238 U 15 16.3 0.7 1.67 3.1 0.195 



Table 2: Comparison of the (3 and 7 band energies (in keV) and E2 transitions (in 



eb) with the sdg-IBM calculations in the actinide region. The data are from |2T| . 



Nucleus 




Ep 




<2 fi \\T(E2)\\0 g > 


<2 7 ||T(£2)||0 9 > 


Cal. 


Exp. 


Cal. 


Exp. 


Cal. Exp. 


Cal. Exp. 


232 Th 


731 


730 


786 


785 


0.27 0.33 ± 0.07 


0.40 0.37 ± 0.07 


234TJ 


804 


809 


872 


927 


0.21 <0.25 


0.35 0.36±0.08 


236 JJ 


903 


919 


978 


958 


0.22 


0.37 


238 JJ 


949 


993 


1040 


1060 


0.22 0.24± 0.05 


0.37 0.37 ± 0.04 
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Figure captions 



Fig. 1. Sytematic study of moment of inertia (a), its spin dependence (c,d), and the 
boson cutoff effect in E2 transitions (d) in the sd-IBM. The quadrupole parameter q 
is normalized to 1 in the SU(3) limit, and the <i-boson energy parameter 77^ = Ed/Nn 
is varied from 0-2 in 10 equal steps. 

Fig. 2. Same as Fig. 1 but in the sdg-TBM.. The quadrupole parameters ^22? Q241 Q44 
are scaled from their SU(3) values with a single factor q, rjd = and the g-boson 
energy parameter r\ g = e g /Nn is varied from 0-10 in 20 equal steps. 

Fig. 3. Comparison of the experimental (circles) and calculated (solid lines) ground 



band energies El/L(L + 1) (in keV) in the actinide nuclei. The data are from [21 . 



Fig. 4. Comparison of the experimental (circles) and calculated (solid lines) yrast 
E2 transitions in the actinide nuclei. The data are from [§, 21]. 
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